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Abstract 

Microtubules are unique cytoskeletal structures that have structural subunits 
of ap tubulin. Taxol is a typical microtubule stabilizing drug. The epothilones 
are other natural products with similar mechanism of action totaxol. Despite 
the highly conserved nature of |3-tubulin, some organism like Saccharomyces 
cerevesia (Sxerevesia) is resistance to taxol, but sensitive to epothilones. In 
order to find differences in sensitivity of yeast tubulin to these molecules, we 
investigated binding mode of the taxol and epothilone A to yeast tubulin 
using molecular modeling. The multiple sequence alignment of p-tubulin of 
different species was performed using ClustalW2. Protein structure of yeast p- 
tubulin was constructed with Swiss Model 8.05 by using 1TVK. Modeled tubulin 
was superimposed with PyMol onUFF for comparison of three-dimensional 
structure of two proteins. Our results showed that one of the most interesting 
differences in binding mode of these molecules is residue 227. The His227 in 
bovine makes a hydrogen bond by means of its 5-nitrogen with epothilone A 
and by means of its s-nitrogen with taxol. The Asn227 of yeast can play role of 
the 5-nitrogen of imidazole ring of H227, but not of e-nitrogen of it. So yeast 
tubulin in contrast to taxol can interact with epothilone A. Due to con- 
servation of essential residues for binding (T274, R282 and 0292), epothilone A 
in comparison with taxol can tolerate the interchange in the binding pocket 
(R276I). Our findings may be of a great aid in the rational design of anti- 
tumor agents that bind to the taxol binding region of tubulin. 
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Introduction 

Microtubules are unique cytoskeletal struc- 
tures that have structural subunits of ap tubu- 
lin, which are important for intracellular 
transport and cell division of all eukaryotes 
An effective and validated target for tumor 
chemotherapeutic agents is the microtubule 
cytoskeleton. These anticancer drugs can be 



divided into two main classes: those that 
inhibit microtubules assembly and those that 
stabilize microtubules. 

Taxol is a typical microtubule stabilizing 
drug of the first group. Other natural products 
with the similar mechanism of action to taxol 
including epothilones, discodermolide, eleu- 
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therobin, sarcadictyin and laulimalides, have 
been isolated and are in various stages of pre- 
clinical/clinical development (2) . 

In experiment with bovine brain tubulin, 
epothilones competitively inhibit the binding 
of taxol to microtubules (3 ' 4) . This data sug- 
gested that there is a common binding site or 
overlapping site for two compounds, despite 
having unrelated structures. In comparison 
with taxol, epothilones are more effective 
than taxol against some multidrug resistance 
cell lines and are more water soluble than 
taxol (3) . 

Despite the highly conserved nature of (3- 
tubulin across the phyla, organisms show dif- 
ferent degrees of susceptibility and resistance 
to the different groups of antimicrotubule 
drugs. For example, many protist and fungal 
tubulins are only slightly affected by taxol (5) . 
However, they are structurally normal and can 
even co-assemble with mammalian tubulins 
(6) . It has been reported that sensitivity to this 
molecule is controlled by relatively few resi- 
dues. One study (7) made a taxol-binding site 
in yeast P-tubulin, which is normally taxol- 
resistant, by replacing five residues with their 
mammalian counterparts (A19K, T23V, G26D, 
N227H, and Y270F). 

According to the high degree of conserva- 
tion among tubulin proteins and the relative 
simplicity of the yeast microtubule cytoskel- 
eton, S.cerevesia is a suitable organism for 
studying the microtubule cytoskeleton. Tubu- 
lin of the budding yeast S.cerevesia shares 
76% sequence identity with bovine tubulin, 
but yeast tubulin isn't sensitive to taxol. In 
contrast, epothilones can promote the assem- 
bly and stabilization of yeast microtubules. So 
the sequence and structure differ in yeast tu- 
bulin and bovine brain tubulin that allow 
strong binding of epothilones, but not of taxol 

(8-11)^ 

In view of the fact that S.cerevesia has only 
one gene (77752) encoding P-tubulin, this or- 
ganism was chosen for our study. In our 
studies, by using molecular modeling we de- 
termined the differences in binding inter- 
actions with yeast P-tubulin among the micro- 



tubule-stabilizing agents, taxol and epothi- 
lones. Understanding of the binding mode 
shared by taxol and other agents with similar 
taxol-like activity (such as epothilones) will 
be a great aid in the rational design of anti- 
tumor agents that bind to the taxol binding 
region of tubulin. 

Materials and Methods 

Multiple sequence alignment 

Protein sequence alignments and compare- 
sons were performed using a Basic Local 
Alignment Search Tool (BLAST) program, 
blastp, against database specification of non- 
redundant proteins which were available at 
the National Center for Biotechnology Infor- 
mation (NCBI) Web server, [http://www.ncbi. 
nih.gov/blast/] (12) . Amino acid sequences 
were obtained from NCBI Web server, [http:// 
www.ncbi.nlm.nih.gov/] . Multiple sequence 
alignments were performed using ClustalW2 
[http://www.ebi.ac.uk/clustalw/], available at 
the European Bioinformatics Institute (EBI) 
Web server (13) . Consensus of amino acid se- 
quence was obtained from Boxshade available 
at the European Molecular Biology Network 
(EMBnet) Web server[http://www.ch.embnet. 
org/]. 

Known 3D structure of proteins 

Three-dimensional protein structures of 
bovine brain tubulin in complex with taxol 
(1JFF) and epothilone A (1TVK) were obtained 
from the Protein Data Bank (PDB) [http// 
www.rcsb.org/pdb/] (1415) . 

Molecular modeling 

Protein structure of yeast P tubulin was 
constructed with Swiss Model 8.05[http:// 
www.expasy.org/swissmodel/] (1619) by using 
a bovine brain tubulin in complex with epo- 
thiloneA (1TVK). The sequence identity be- 
tween yeast P tubulin and bovine brain tubu- 
lin are 76%. The yeast tubulin structure model 
was quite similar to the bovine brain tubulin 
crystal structure. The Root Mean Square 
Deviation (RMSD) for 1704 polypeptide back- 
bone atoms was 0.11 A as calculated with 
Swiss-Pdb Viewer version 4.0.1[http://www. 
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expasy.org/spdbv/] (20) . We used QMEAN Z- 
score as an estimate of the model quality. The 
reliability of the modeled structure was 
checked by PROCHECK software online. 

Molecular modeling visualization and analysis 

Modeled yeast tubulin was superimposed 
with PyMol on the bovine (3-tubulin in com- 
plex with taxol (1JFF) for comparison of 
three-dimensional structure of two proteins. 
Structural presentation of 1JFF and 1TVK and 
modeled yeast tubulin was made by using 
Swiss-Pdb Viewer and PyMOL (Version 0.99) 
[http://pymol.sourceforge.net] programs (21) . 
Cartoon drawings of the structures were ob- 
tained using PyMOL (21) . Studies and compare- 
sons of taxol and epothiloneA binding inter- 
action with yeast tubulin were performed by 
using programs, Swiss-Pdb Viewer and PyMOL. 

Results 

The multiple sequence alignment of (3-tubu- 
lin of different species is shown in figure 1. 
Consensus sequences obtained from the Box- 
shade program is shown. According to this 
alignment (3-tubulin is highly conserved in 
different eukaryote. 

Tubulin has binding site for GTP that high- 
ly conserved among the species. Tubulin like 
other nucleotide-binding proteins belongs to a 
class of regulatory proteins called G proteins. 
Consensus sequence GQC, KGHYTEG, DNE 
AL and glycine rich sequences SLGGGTGSG 
MG are binding site for phosphoryl group of 
nucleotide and are observed in other nucleo- 
tide binding proteins (Figure 1) (22) . Residues 
in the 15-27, 212-231, 270-285 and 358-372 
region of (3-tubulin form hydrophobic sites for 
taxol (23) and residues in the 212-231 and 270- 
282 regions make contact with epothilone (24) . 

In these regions, we observed some differ- 
ences between yeast (3-tubulin and (3-tubulin 
of other species. For example, the yeast tubu- 
lin instead of K19, V23, D26, H227, A231, 
F270, R276 and R359 has A, T, G, N, S, Y, I 
and Q respectively (Figure 1). 

The modeled yeast (3-tubulin superimpos- 
ing on the bovine (3-tubulin with a RMSD of 
0.1 1 A and E-value of 4.2e-45 showed that the 
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Figure 1. Multiple alignment of P-tubulin of different spe- 
cies: a) Oryza sativa, b) Plasmodium falciparum, c) Salmo- 
salar, d) Bostaurus, e) Drosophila melanogaster, f) Homo 
sapiens, g) Cyathostomumcatinatum, h) Gallus gallus and i) 
Saccharomyces cerevesia. Key residues for taxol binding 
(red), epothilone A binding (orange) and common residues 
for binding of taxol and epothilone A (green) are highlighted 

overall structure is similar (Figure 2). The 
QMEAN Z-score directly shows how many 
standard deviations the model's QMEAN score 
differs from expected values for native or ex- 
perimental structures {e.g. unexpected solvent 
accessibility, backbone geometry, interatomic 
packing, etc.) (25) .Good models reach QMEAN 
Z-scores comparable to experimental struc- 
tures (0<mean Z-score<2). 

According to QMEAN Z-score analysis, 
modeled structure receives a Z-score of -1.8. 
PROCHECK was used to analyze the Rama- 
chandran plot of the model by comparing it 
with the Ramachandran plots of 163 represen- 
tative structures solved at high resolution of at 
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M loop 

Figure 2. Superimposition of backbone atoms of modeled 
yeast P-tubulin (light magenta) on bovine P-tubulin 
(cyan). Arrows indicates the position of the M loop and 
residue 276 in the P-tubulin 

least 2.0 A. The Ramachandran plot analysis 
of the model structure showed that 97.20% of 
the residues are in favor of additional allowed 
regions. The model was thus evaluated to be 
of good quality based on these results (26) . The 
most striking difference lies in the M loop, 
which is the taxol binding site. One of the 
most important residues in M loop that affects 



binding of taxol to tubulin is residue 276. 
Comparison of interaction between taxol or 
epothilone and residue 276 of yeast and 
mammalian P-tubulin is presented in figure 3. 

Analysis of binding interaction of bovine (3- 
tubulin and taxol reveals that 17 residues from 
the hydrophobic pocket of taxol binding site 
made direct contact with toxol molecule in- 
cluding V23, D26, L215, L217, H227, L228, 
A231, S234, F270, P272, L273, T274, S275, 
R276, P358, R359, G360 (Figure 4A). Seven of 
these residues are different in yeast P-tubulin 
(Val23Thr, Asp26Gly, His227Asn, Ala231Ser, 
Phe270Tyr, Arg276Ile, Arg359Gln replacement) 
(Figure 4B). 

Binding mode analysis of bovine P-tubulin 
and epothilone A showed that 6 residue from 
the epothilone A binding site made direct con- 
tact with epothilon A molecule including 
H227, A231, T274, R276, R282,Q292 (Figure 
4C). Yeast P-tubulin instead of H227, A231 
and R276 has N, S and I respectively (Figure 
4D). 

Molecular modeling shows that another 
important residue for binding of taxol and 
epothilone is position 227. Mode of interact- 




Figure 3. A) Residue 276 of bovine tubulin and its interaction with taxol, B) Residue 276 of yeast tubulin and its interaction 
with taxol, C) Residue 276 of yeast tubulin and its interaction with epothilone A and D) Residue 276 of bovine tubulin and its 
interaction with epothilone A 
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Figure 4. Comparison of the binding mode of paclitaxel and epothilones in yeast and bovine tubulin. Ligands (gray), binding-site 
key residues (greencyan) and non-key residues (slate) for bovine tubuline, and), binding-site key residues (salmon) and non-key 
residues (slate) for yeast tubuline that e are shown in stick representation. Hydrogen bonds are presented as yellow long dashes 
and hydrophobic interactions are shown as red square dots. The remaining part of the P-tubulin secondary structure is rendered 
with cartoon presentation using PyMOL. A) Taxol as found in the electron crystallography model (PDB: 1JFF). B) Taxol that 
interacts with modeled yeast tubulin. C) Epothilone A as found in the electron crystallography model (1TVK). D) Epothilone A 
that interacts with modeled yeast tubulin 



tion of taxol and epothilone with this residue 
is compared in figure 5. This comparison 
shows that these molecules have different 
mode of interaction with this residue. 

Discussion 

Our alignment showed the high identity 
and homology in sequence of P-tubulin of 
different species. We didn't see any gap 
except in C terminal. At C-terminal of P-tubu- 
lin, we observed extensive molecular hetero- 
geneity. These C terminal regions are highly 
variable among the isotype within a species 
and among species. This region is important 
for anchoring Microtubule Associated Protein 
(MAP) in place (27) . The diversity of C- 



terminal may propose the mechanism for 
isotype-specific MAP binding. 

Molecular modeling shows that M loop of 
yeast P-tubulin in comparison with bovine P- 
tubulin locates far away taxol molecule. M 
loop is the main component in contacts be- 
tween protofilaments, interacting with struc- 
tures in the H1-S2 and H2-S3 loops of the 
adjacent monomer (28) . S.cerevesia exhibit 
some nonconservative substitutions in the M- 
loop compared to other species (Figure 1). 
The interesting substitution is at position 276. 
The residue R276 in the M loop of bovine P- 
tubulin, which is a positively charged residue, 
is replaced by an isoleusine in yeast. This 
position exists near a predicted bend in the M- 
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Figure 5. Possibility of hydrogen bonding between residue 227 of tubulin [Histidine in bovine (cyan) and Asparagine in 
yeast (magenta)] and A) oxygen atom of 3 '-benzamidophenyl group of taxol, B) nitrogen atom of thiazol ring of 
epothilone A 



loop, and an isoleusine residue here may 
change the shape of M-loop and its interaction 
with taxol. Because of the flexible long side 
chain and its surface location, arginine can 
interact with taxol. 

There are many oxygen groups in taxol 
molecule that can interact and make hydrogen 
bond with positively charged goanidino group 
of arginine (Figure 3A). This interaction can 
stabilize and keep the molecule of taxol in its 
hydrophobic pocket. In contrast, neutral iso- 
leusine cannot make such hydrogen bond and 
thus cannot stabilize taxol in its pocket 
binding (Figure 3B). The binding of epothil- 
one A do not affect by this change so the 
yeast P-tubulin is yet sensitive to this drug. 
The epothiloneA molecule can interact with 
yeast tubulin because of existence of im- 
portant hydrogen bonds between epothilone A 
and yeast tubulin (Figures 3C and 3D). 

The isopropyl group of V23 is near the 
C3 '-benzamidophenyl ring in the bovine P- 
tubulin (Figure 4A) (29) . Yeast P-tubulin has a 
theronine residue instead of valine. Substi- 
tution of non-polar valin by polar threonine 
leads to distort the positioning of the hydro- 
phobic phenyl ring (Figure 4B). In the bovine 
P-tubulin, the methylene group of K19, E22, 
D26 and E27 interact with the C3 '-benz- 
amidophenyl ring and stabilize it (28) . In yeast 
P-tubulin, replacement of K19 and D26 by 
alanin and glycin respectively reduce the 
methylene groups that interact with the C3- 
benzamidophenyl ring. The size of epothilone 



A molecule is smaller than taxol molecule so 
epothilone A only fill half of taxol binding 
pocket in P-tubulin. Thus, residues in 15-27 
region of P-tubulin including K19, V23 and 
D26 which interact with taxol, do not play an 
important role in epothilone binding (Figure 
4C) (29) . We measured distance between hy- 
droxyl group of T23 of yeast P tubulin 
predicted model and the C3 '-benzamido- 
phenyl ring of taxol molecule, and distance 
between OH group of T23 of yeast modeled P 
tubulin and thiazol ring of epothilone. These 
distances were 5.7 A for taxol and 9.88 A for 
epothilone A. So the polarity of hydroxyl 
group of threonine can affect the binding of 
taxol, but not of epothilone A. 

In a mutagenesis study, Himes et al re- 
ported that epothilones despite taxol can 
tolerate the change made at position 227 and 
bind equally well to wide-type (N227) and the 
mutated (H227) yeast tubulin (7) . 

According to the model, the imidazol ring 
of H227 of bovine P-tubulin is situated be- 
tween C2-benzoyl phenyl and C3 '-benz- 
amidophenyl ring of taxol molecule (Figure 
4A). The s-nitrogen of imidzol ring of his- 
tidine 227 can make hydrogen bond with 
oxygen atom of C3 ' benzamidophenyl ring of 
taxol, and thus stabilizes this drug in its 
hydrophobic pocket (Figure 5 A). The thiazol 
ring of epothilone A and the imidazol ring of 
H227 of bovine P-tubulin are located in a 
face-to-face orientation (Figure 5B) (30) . The 
5-nitrogen of imidazole ring of H227 can 
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make hydrogen bond with nitrogen atom of 
thiazol ring of epothilone A (Figure 5B). 

In our yeast (3-tubulin model, asparagine 
was situated position 227 (Figure 5A) and 
because of its distance, 5-nitrogen atom of 
N227 cannot play the role of 8-nitrogen of 
imidzol ring of H227, thus eliminates the 
hydrogen bonds (Figure 5 A). In contrast, the 
5-nitrogen atom of N227 in yeast tubulin can 
mimic the role of 5-nitrogen atom of histidine 
and makes hydrogen bond with thiazol ring of 
epothilone A (Figure 5B). 

According to this model, the residue F270 
in bovine (3-tubulin was located in a hydro- 
phobic basin (Figure 4A) (28) . This residue can 
make van der Waals contact with C3 '-phenyl 
ring and C4- acetyl group of taxol molecule. 
The epoxide ring at position CI 2, CI 3 of 
epothilone A accommodates the C3 '-phenyl 
ring and C4- acetyl group of taxol (31) , and is 
located at the entrance of the hydrophobic 
basin, but in contrast totaxol, it doesn't enter 
this pocket. The measured distance between 
the C3 '-phenyl ring of taxol and F270 was 3.8 
A. This distance for epoxide ring of epothi- 
lone A and same residue was 8.96 A. In our 
yeast (3-tubulin model, tyrosine and serine 
were located at positions 231 and 270. These 
hydrophilic side chains distort the van der 
Waals contacts of hydrophobic basin with 
C3 '-phenyl ring and C4- acetyl group of taxol 
molecule. The residue F270 is a key residue 
for taxol binding, but not to the epothilone A. 
So the change of this residue in contrast to 
taxol does not affect the binding of epothilone 
to yeast (3-tubulin. The cytotoxicity data show 
that the Phe270Val mutation has only a 3 -fold 
effect on the sensitivity of the epothilone A 
compared with a 27-fold change for the taxol 
(32) , confirm this finding. 

The position 359 in bovine (3-tubulin is 
occupied by arginine that makes hydrogen 
bond with taxol molecule (Figure 4A). In 
yeast (3-tubulin glutamine locates at position 
359 and distorts this hydrogen bonding. In 
agreement with our finding, Henriquez et al 
reported that sequence replacement of Arg 
359Ala in Acanthamoeba tubulin lead to 



resistance to taxol (33) . The distances of this 
residue with taxol and epothilone A was 
measured 3.2 and 8.4^, respectively. Accord- 
ing to this finding, epothilone is so far to 
make direct hydrogen bound with this residue 
so that its binding is not affected by changing 
this residue. 

According to the model, the 7-OH group of 
epothilone and oxetane ring of taxol made 
hydrogen bonding with T274. The residue 
T274 is important for binding of both of taxol 
and epothilone A and it is confirmed by muta- 
genesis study which reported that a Thr274Ile 
replacement had a determined effect on 
activity (34) . This effect is more emphasized on 
epothilone than on taxol activity. This residue 
is conserved in yeast so the epothilone A can 
made one of the most important interaction 
with its (3-tubulin. 

The model also showed that residue argi- 
nine at position 282 is essential for binding of 
epothilone (Figure 4B), but not of taxol. 
Again, this residue is conserved in yeast 
thereby epothilone A can bind to its (3-tubulin. 
This positively charged residue made electro- 
static contacts with the 7-OH group of epothi- 
lone A. On the other hand, R282 can form hy- 
drogen binding with T274 and thereby stabil- 
izing conformation of M loop for interaction 
with epothilone A. In contrast, taxol doesn't 
make such network of electrostatic inter- 
actions and this residue isn't a key residue for 
its binding and again mutagenesis study 
agrees with this fact (24) . 

The other residue that is important for bind- 
ing of epothilone, but not to the taxol is Q292 
and again this residue is conserved in yeast, 
so this organism is sensitive to epothilone A 
(Figure 4B). This residue locates near the M 
loop from and forms hydrogen binding with 
Leu273 adjacent to T274 (15) . Mutagenesis 
studies showed that Gln292Glu replacement 
distorts this hydrogen bonds network of M 
loop so inhibits epothilone binding (34) . 

Conclusion 

Analysis of S.cerevesia (3-tubulin has 
shown that failure in interaction with the C3- 
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benzamidophenyl and phenyl rings and stabil- 
izing them, because of Val23Thr, Asp26Gly, 
His227Asn, Ala231Ser, Phe270Tyr replace- 
ments, is the main reason for resistance to 
taxol. Study of modeled yeast tubulin re- 
vealed that ability of epothilone to stabilize 
yeast tubulin is related to conservation of key 
residues for epothilone binding (including 
T274, R282 and Q292) in this organism. On 
the other hand, epothilone A can tolerate 
some interchange in taxol binding pocket resi- 
dues due to flexibility of epothilone A in 
comparison with taxol. 

Our direct comparison of epothilone bind- 
ing and taxol binding in yeast tubulin can sup- 
port development of novel natural and syn- 
thetic microtubule stabilizing agents with 
higher affinity for tubulin and lower sensitiv- 
ity toward development of resistance of can- 
cerous cells. This would lead to make anti- 
cancer drugs with an improved efficiency and 
fewer side effects. 

Acknowledgement 

Special thanks are extended to Dr. Ahmad 
Reza Shahverdi and Dr. Mojtaba Panjehpour 
for their helpful suggestions. This research 
was supported by the Bioinformatics Re- 
search Center of Isfahan University of Med- 
ical Sciences. 

References 

1. Hyams JS, Lloyd CW. Microtubules. In: Harford 
JB (eds). Modern Cell Biology. Vol. 13. New 
York:Wiley-Liss;1994. 

2. He L, Orr GA, Horwitz SB. Novel molecules that 
interact with microtubules and have functional ac- 
tivity similar to Taxol. Drug Discov Today 200 1;6 
(22):1 153-1 164. 

3. Bollag, DM, McQueney PA, Zhu J, Hensens O, 
Koupal L, Liesch J, et al. Epothilones, a new class 
of microtubule- stabilizing agents with a taxol-like 
mechanism of action. Cancer Res 1995;55(11): 
2325-2333. 

4. Kowalski R J, Giannakakou P, Hamel E. Activities 
of the microtubule-stabilizing agents epothilones A 
and B with purified tubulin and in cells resistant to 
paclitaxel (Taxol®). J Biol Chem 1997;272(4): 
2534-2541. 



5. Gull K. Protist tubulins: new arrivals, evolutionary 
relationships and insights into cytoskeletal func- 
tion. Curr Opin Microbiol 2001;4(4):427-432. 

6. McKean PG, Vaughan S, Gull K. The extended 
tubulin superfamily. J Cell Sci 2001;114(Pt 15): 
2723-2733. 

7. Gupta ML, Bode CJ, Georg GI, Himes RH. 
Understanding tubulin-taxol interactions: mutations 
that impart taxol bining to yeast tubulin. Proc Natl 
Acad Sci USA 2003;100(1 1):6394-6397. 

8. Little M, Seehaus T. Comparative analysis of tubu- 
lin sequences. Comp Biochem Physiol B 1988;90 
(4):655-670. 

9. Burns RG. Alpha-, beta-, and gamma- tubulins: se- 
quence comparisons and structural constraints. Cell 
Motil Cytoskeleton 1991;20(3): 181-189. 

10. Brans G, Louie KA, Botstein D. Yeast proteins as- 
sociated with microtubules in vitro and in vivo. 
Mol Biol Cell 1992;3(l):29-47. 

11. Bode CJ, Gupta ML, Reiff EA, Suprenant KA, 
Georg GI, Himes RH. Epothilone and paclitaxel: 
unexpected differences in promoting the assembly 
and stabilization of yeast microtubules. Biochem- 
istry 2002;41(12):3870-3874. 

12. Altschul SF, Madden TL, Schaffer AA, Zhang J, 
Zhang Z, Miller W, et al. Gapped BLAST and PSI- 
BLAST: a new generation of protein database 
search programs. Nucleic Acids Res 1997;25(17): 
3398-3402. 

13. Thompson JD, Higgins DG, Gibson TJ. CLUSTAL 
W: improving the sensitivity of progressive mul- 
tiple sequence alignment through sequence weight- 
ing, position-specific gap penalties and weight 
matrix choice. Nucleic Acids Res 1994;22(22): 
4673-4680. 

14. Lowe J, Li H, Downing KH, Nogales E. Refined 
structure of alpha beta-Tubulin at 3.5 A resolution. 
J Mol Biol 2001 ;3 13(5): 1045-1057. 

15. Nettles JH, Li H, Cornett B, Krahn JM, Snyder JP, 
Downing KH. The binding mode of epothilone A 
on alpha, beta-tubulin by electron crystallography. 
Science 2004;305(5685):866-869. 

16. Arnold K, Bordoli L, Kopp J, Schwede T. The 
SWISS-MODEL Workspace: A web-based envir- 
onment for protein structure homology modelling. 
Bioinformatics 2006;22(2): 195-201. 

17. Kiefer F, Arnold K, Kiinzli M, Bordoli L, Schwede 
T. The SWISS-MODEL Repository and associated 
resources. Nucleic Acids Res 2009;37(Database 
issue):D387-D392. 

18. Schwede T, Kopp J, Guex N, Peitsch MC. SWISS- 



Avicenna Journal of Medical Biotechnology, Vol. 3, No. 4, October-December 2011 



Abbari V, et al 



MODEL: an automated protein homology-model- 
ing server. Nucleic Acids Res 2003;31(13):3381- 
3385. 

19. Guex N, Peitsch MC. SWISS-MODEL and the 
Swiss-Pdb viewer: an environment for comparative 
protein modelling. Electrophoresis 1997; 18(1 5): 
2714-2723. 

20. Kabsch W. A discussion of the solution for the best 
rotation to relate two sets of vectors. Acta Cryst 
1978;34(Part 5):827-828. 

21. Delano WL. The PyMOL molecular graphics sys- 
tem. DeLano Scientific; San Carlos, CA, USA: 
2005. 

22. Sternlicht H, Yaffe MB, Farr GW. A model of the 
nucleotide-binding site in tubulin. FEBS Lett 1987; 
214(2):226-235. 

23. Xiao H, Verdier-Pinard P, Fernandez-Fuentes N, 
Burd B, Angeletti R, Fiser A, et al. Insights into the 
mechanism of microtubule stabilization by Taxol. 
Proc Natl Acad Sci USA 2006;103(27):10166- 
10173. 

24. Giannakaou P, Gussio R, Nogales E, Downing KH, 
Zaharevitz D, Bollbuck B, et al. A common phar- 
macophore for epothilone and taxanes: Molecular 
basis for drug resistance conferred by tubulin mu- 
tations in human cancer cells. Proc Natl Acad Sci 
USA 2000;97(6):2904-2909. 

25. Benkert P, Biasini M, Schwede T. Toward the esti- 
mation of the absolute quality of individual protein 
structure models. Bioinformatics 201 1;27(3):343- 
350. 

26. Laskowski RA, MacArthur MW, Moss DS, Thorn- 
ton JM. PROCHECK: a program to check the ste- 
reochemical quality of protein structures. J Appl 
Cryst 1993;26(2):283-291. 



27. Amos LA. Focusing-in on microtubules. Curr Opin 
Struct Biol 2000;10(2):236-241. 

28. Nogales E, Whittaker M, Milligan RA, Downing 
KH. High-resolution model of the microtubule. 
Cell 1999;96(l):79-88. 

29. Synder JP, Nettles JH, Cotnett B, Downing KH, 
Nogales E. The binding conformation of Taxol in 
(3-tubulin: A model based on electron crystallo- 
graphy density. Proc Natl Acad Sci USA 2001;98 
(9):5312-5316. 

30. Reese M, Sunchez-Pedregal VM, Kubicek K, Mei- 
ler J, Blommers MJJ, Griesinger C, et al. Structural 
basis of the activity of the microtubule-stabilizing 
agent epothilone a studied by NMR spectroscopy in 
solution. Angew Chem Int Ed 2007;46(1 1):1864- 
1868. 

31. Forli S, Manetti F, Altmann KH, Botta M. Evalu- 
ation of novel epothilone analogues by means of a 
common pharmacophore and a QSAR pseudore- 
ceptor model for taxanes and epothilones. Chem 
MedChem 2010;5(l):35-40. 

32. Giannakakou P, Sackett DL, Kang YK, Zhan Z, 
Buters JT, Fojo T, et al. Paclitaxel-resistant human 
ovarian cancer cells have mutant beta-tubulins that 
exhibit impaired paclitaxel-driven tubulin polymer- 
ization. J Biol Chem 1997;272(27):171 18-17125. 

33. Henriquez FL, Ingram PR, Muench SP, Rice DW, 
Roberts CW. Molecular basis for resistance of 
Acanthamoeba tubulins to all major classes of anti- 
tubulin compounds. Antimicrob Agents Chemother 
2008;52(3):1 133-1 135. 

34. He L, Yang CPH, Horwitz SB. Mutations in (3- 
tubulin map to domains involved in regulation of 
microtubule stability in epothilone-resistant cell 
lines. Mol Cancer Ther 200 1 ; 1 ( 1 ) : 3 - 1 0. 



Avicenna Journal of Medical Biotechnology, Vol. 3, No. 4, October-December 2011 



